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Several thermodynamic quantities within the planar t — J model are calculated using the T > 
Lanczos method on clusters of up to 26 sites. Hole density Ch(fJ-,T) shows a non- Fermi liquid 
behavior as a function of T and suggests a transition from small to large Fermi surface at cu ~ 0.15. 
Specific heat reveals a maximum at the exchange energy-scale up to = 0.2, where it becomes 
almost T- independent for T > 0.15 t. At constant T the entropy has maximum for Ch ~ 0.15, with 
large values at low T, consistent with experiments on cuprates. In the underdoped regime the spin 
susceptibility Xo{T) exhibits a maximum at finite T — T* , with T* decreasing with doping, and 
disappearing for Ch > 0.15. 



PACS numbers: 71.27.+a, 75.40. Cx, 71.10.+X 

The normal state of the superconducting cuprates, in 
particular the character of the low-energy excitations, 
still lacks proper understanding. Although there ap- 
pears to be a consolidation of experimental results prob- 
ing their transport 0] , magnetic [0 , thermodynamic || , 
and single-electron properties, on the theoretical side 
several important questions about the nature of the cor- 
responding quantum liquid remain open. While un- 
doped cuprates are well understood in terms of a two- 
dimensional (2D) 5=1/2 antifcrromagnct (AFM), it is 
less clear what changes are induced on doping the AFM 
with mobile holes. Besides the questionable existence of 
quasiparticles (QP) in the sense of some type of Fermi 
liquid, the character of the QP and the associated Fermi 
surface (FS) in doped AFM are also controversial. While 
at low doping theoretical studies |^|j6) as well as exper- 
imental findings [Q seem to favor the interpretation in 
terms of hole QP in a rigid band and related hole-pocket 
FS, situation at larger doping is more consistent with 
the usual large electronic FS |,|,f§. Nevertheless, sev- 
eral normal-state dynamical properties of cuprates have 
recently been reproduced by using unbiased numerical 
techniques with strong correlations taken fully into ac- 
count |l(]Jll| • It is the aim of this paper to present results 
on relevant thermodynamic properties of correlated elec- 
trons, which could help in finding a proper concept for 
such systems, but also in establishing the relevant model 
for cuprates. 

In the following we address some of the questions above 
by studying finite-T static properties of the 2D t — J 
model, as a prototype model for strongly correlated elec- 
trons and cuprates in particular 
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Here c| , Ci S are projected fermionic operators, excluding 
the states with doubly-occupied sites. To investigate the 
regime of cuprates we set J/t — 0.3. 



Although there have been so far extensive numeri- 
cal studies of the ground-state properties of the t — J 
model its finite-T properties are less explored. AT 
T > a fruitful technique in addressing several thermo- 
dynamic properties proved to be the high-T series ex- 
pansion [fnjJlSp lj , In the present study we aim at the 
calculation of the finite-T static properties of the model 
Eq. (p]) with a recently introduced method, combining 
Lanczos diagonalization technique and random sampling 
|]i"of , which has been so far applied to the study of several 
dynamical responses |ll]] . We concentrate on calculating 
the quantities, which are related to the expectation values 
of the conserved operators, commuting with the Hamil- 
tonian, Eq.(l). Such operators are e.g. total number of 
electrons N e , total energy H, and S z , the z-component 
of the total spin operator. In terms of them some in- 
teresting thermodynamic quantities can be expressed, in 
particular chemical potential //, entropy S, specific heat 
Cy and the uniform susceptibility xo- The reason to 
limit our treatment only to commuting operators is tech- 
nical, since in this case we can avoid time and memory 
consuming computation of wavefunctions and their scalar 
products. This in turn enables us to study larger clusters 
with up to 20 sites at all doping concentrations. 

We follow the method for the evaluation of T > 
static quantities Jl(]]. To calculate the quantity of the 
type Tr[f{N e ,S z ,H)e-^ H -i J - N ^}, the trace is expanded 
in terms of an arbitrary orthonormal set of basis func- 
tions |n), n = 1 . . .K, with good quantum numbers N e 
and S z - Each state \n) is used as initial function |^q) for 
the M-step Lanczos procedure fL3|] , yielding a subspace 
spanned by functions j = . . . M. The Hamiltonian 
is diagonalized in this subspace to obtain the eigenvalues 
E™ and corresponding eigenvectors ), which are used 
to evaluate (/) 

K M-l 
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where Q = Tr[e~P( H ~ IJ ' N '^] is the grand-canonical par- 
tition function. Clearly, the expression Eq.(||) is exact 
for K = M = N st , where N st is the dimension of the 
complete basis. In practice K,M <C N st are used, with 
states \n) sampled randomly. More detailed discussion of 
the method and tests are given elsewhere In the 

present application it is essential that scalar products 
(n\ipj) are trivially evaluated, since this is the compo- 
nent of the function |0q) in the wavefunction i.e. 
can be obtained from the eigenvectors of the M x M 
tridiagonal matrix without storing </>"}■ Typically also 
small M < 100 is used, thus reorthogonalization of the 
Lanczos functions can be avoided. The computational 
effort in the present case is then equal to a ground-state 
Lanczos procedure, repeated K times. We employ in the 
following typically K ~ 200 — 1000 in every N e sector. 

We have calculated a few static quantities above using 
the approximation Eq. (0) on clusters with N — 16, 18, 
and 20 sites, thus having control over the finite-size ef- 
fects through comparison of the results. We also treated 
the undoped system with N = 26 sites. In general the 
results are quite free of finite-size effects for temperatures 
T > T is (N, c h ), which amounts to T fs ~ 0.1 1 for N = 20 
and Ch < 0.3. 

We first analyse the hole density Ch = 1 — (N e )/N as 
a function of T and /i. In Fig. [l](a) we present curves 
Ch(T) for several fih, where the chemical potential for 
holes is Hh = — /i. Clearly, for /i^ < /i" the ground state 
of the system contains no holes. Here /i° = —1.99 t JlS] ] 
is related to the minimum energy of a single hole added 
into the undoped system. The corresponding curves for 
[ih < fJ'l are denoted with dashed lines. In Fig. [l](a) we 
in general do not find a T 2 dependence of Ch at low T, 
as expected for a normal Fermi liquid, except within the 
extreme overdoped regime cu > 0.3. In particular, in a 
broad range 0.15 < Ch < 0.3 very unusual linear variation 
of Ch(T) is observed at lowest T. Moreover, a remarkable 
feature is a non-monotonous c^(T) dependence as the 
value c/,(0) is varied. There exists a doping concentration 
c* h such that derivative dch/dT at low T is positive for 
curves with Ch(T ~ 0) < c* h and negative for curves with 
Ch(T ~ 0) > c* h . The marginal doping c* h ~ 0.15 seems 
to be system independent, as checked quantitatively for 
systems with N = 16, 18, 20. 

In Fig. 0(b) we plot the variation of c/j with at 
several chosen temperatures T. We include also T = 
0.05 t < Tf s , since in this case it shows less pronounced 
finite-size effects. As T — > clearly Ch(^h < ~~ * 0> 
while Ch{^h > Mft) remains finite. The slope dch/d^ih, 
being proportional to the compressibility n of the hole 
fluid, is finite for T > 0, indicating the absence of the 
phase separation in the system at chosen J/t pa]. On 
the other hand, Fig. |l|(b) reveals that for T — > K is 
increasing and becoming large on approaching ~ 
consistent with large (polaron) mass enhancement of in- 
dividual holes at vanishing doping. Moreover we cannot 



exclude a singular behavior (divergence) of deduced 
from analogous results for the Hubbard model jl6| . 

Next we find remarkable a feature, appearing at the 
hole-density of ~ C/*, where \Xh is essentially pinned 
at the value fih ~ — 1.8 t, i.e. does not vary with 
T. This pinning is active in a wide range of T, and 
is again almost not dependent on the system size. It 
is tempting to interpret the marginal concentration c* h 
as a change of the character of the FS. To establish 
the relation, we still have to rely on arguments which 
apply to the gas of non-interacting fermions. Denot- 
ing by g(e) the density of the single-fermion states, a 
simple Sommerfeld expansion yields that the number 
of particles at fixed chemical potential /i is given by 
c e (T) = c e (T = 0) + (tt 2 /6)(k B T) 2 g'(n). Therefore, the 
prefactor of the T 2 term is just proportional to the deriva- 
tive of the density of states. Indirectly this gives an in- 
formation about the character of the FS, since one would 
plausibly associate g'(fi) > for c e < 1 with a large 
electron FS, and oppositely g'(fi) < with hole-like FS 
or small hole pockets vanishing for c e — > 1. For 2D it 
is straightforward to represent g'{n) more explicitly as a 
line integral 
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over the curve L M = {k; e(k) = /i}. Here 1/K is 
the curvature of (positive for electron-like orbits and 
vice versa), K the unit vector perpendicular to L M , and 
m 1 = d 2 e(k)/dkdk the inverse effective mass tensor. 
Apart from the effective mass term, this is similar to the 
expression for the Hall resistivity in terms of the cur- 
vature of the FS @. From Eq. (§) it follows that at 
least in the region of the /c-space, where the effective- 
mass tensor is positive-definite, g'(p) > implies also 
that the average curvature K^ 1 of the FS L e is positive. 
Although the observed non-quadratic T-dependence in 
Fig. Q(a) questions above arguments, we may still inter- 
pret dch/dT < (i.e. dc e /dT > 0) in Fig. ||(a) as an 
indication for g'(/i) > 0, i.e. positive average curvature 



of the FS for c h > 



This in turn implies a transi- 



tion from the hole-pocket picture [^-Q, existing at low 
doping, to an electron- like FS [|]|J at Ch ~ c* h - 
Next we consider the entropy 



S = k B lnQ+((H}-fi(N e })/T. 



(4) 



As a test we have checked and found agreement of our 
entropy data with available results from the high-T sc- 
ries |l^]. In Fig. U we show the doping dependence of 
the entropy density (per site) at different T < J and 
N = 16 — 20. As expected, finite size effects are most 
pronounced at lowest T = 0.1 t ~ Tf s . Nevertheless, at 
all T < J the entropy has a broad maximum at Ch ~ 0.15, 
indicating the highest density of many-body states in this 
'optimum' doping regime. 
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The magnitude of the entropy in the 'optimum' regime 
appears very large, i.e. at T = 0.1 t = J/3 the en- 
tropy per site is s ~ 0.39 ks, which is almost 40% of 
s(T = oc) for the same doping Ch, although T <C J, t. In 
other words, introducing the degeneracy temperature by 
s(Tdeg) = s(T = oo)/2, for intermediate dopings of Ch ~ 
0.15 we get Tkeg ~ 0.17 t, being very small in comparison 
with any resonable QP bandwidth W, e.g. W = 8 t of 
the non-interacting electrons. It is remarkable that the 
entropy of such magnitude has been deduced from the 
electronic specific heat measurements in YP^CuaOg+x 
materials ||]. E.g., for the optimally doped material with 
x = 0.97 at T = 300 K it was found As = 0.35 k B per pla- 
nar copper site (As = 0.7 k B per unit cell), relative to the 
undoped x = sample. We find the corresponding value 
s(c h = 0.15) - s(c h = 0) ~ 0.30 k B at T = 0.1 t ~ 450 K 
(assuming t = 0.4 eV @). 

In Fig. H we present the T-dependence of the specific 
heat C v = T(dS/dT)^/N at different dopings. For the 
undoped AFM we have not observed any appreciable size 
dependence on systems with 16 - 26 sites, and our results 
seem to be even superior to those obtained by other meth- 
ods |llj. Cy is strongly T-dependent in the observed 
interval T = — 3.3 J, with a maximum at T ~ 2 J/3, 
and approximately Cy oc T 2 at low-T. As the system is 
doped, Cy(T) still exhibits a maximum, which is how- 
ever strongly suppressed and gradually moves to lower 
T with increasing Ch . The peak can be attributed to the 
thermal activation of the spin degrees of freedom with an 
exchange energy scale persisting in the doped system, as 
observed already in dynamical spin correlations p"l| ] . The 
presence of the exchange energy scale however disappears 
in the 'overdoped' regime Ch > 0.3. 

It is characteristic (and consistent with the vanishing 
role of J) that in the 'optimally' doped regime Ch ~ 0.2 
we find CV(T) ~ const for 0.15 t < T < t, being far 
from the Fermi liquid behavior. The anomalous behav- 
ior in this regime implies Q oc T, as follows from Eq. (Q), 
as well as a constant (fiat) density of many-body states 
D(E) ~ const, as defined through fi = / D(E) e - l3E dE, 
at least in a broad intermediate region of energies. It 
seems natural to connect this phenomenon with the ob- 
served universality of dynamical response in the optimum 
doping regime jy] . 

Only at extreme overdoping Ch > 0.5 (low electron 
density) we were able to see indications of the Fermi- 
liquid-like Cv oc T behavior in a broad interval T < t. 
On the other hand, we have indications for a different 
Cv ~ 7T behavior at lowest Tf s < T < 0.15 t, as de- 
duced from a drop of Cv in Fig. [jj, which would be closer 
to the experimental results || and the QP picture. Such 
interpretation would be again difficult, since it would re- 
quire Tdog ~ 800 K. On the other hand, it is puzzling 
that obtained 7 in the low T < 0.15 t regime can be 
matched with the value for non-interacting fermions at 
c e < 1 in a tight-binding band renormalized only by a 



modest factor t/t ~ 2.5, as noted previously ||. 

As last we consider the uniform spin susceptibility 
Xo = {(S z ) 2 )/Nk B T, which is shown in Fig. || at sev- 
eral Ch- Again for = our results agree with other 
calculations E3] down to T ~ 0.3 J. In the undoped 
system \o exhibits a maximum at T = T* ~ J, which 
gradually shifts to lower T with doping and finally disap- 
pears at Ch > 0.15. xo at T — > increases with doping, in 
agreement with experiments Q and some previous results 
|l4],[l3]]. The existence of such maximum at T — T* in 
undoped cuprates has been interpreted in terms of a pseu- 
dogap At doping Ch = 0.15 we observe a monotonous, 
Pauli-like xo f° r T < 0.2 t, which could signify an onset 
of a QP low-T behavior. On the scale of temperatures 
of the order of t, however, we don't find the Fermi liquid 
behavior xa — const up to Ch > 0.6. 

To conclude, we have calculated several thermody- 
namic properties of strongly correlated electrons within 
the t — J model and found essentially non-Fermi liquid 
behavior in a wide range of doping concentrations, up 
to Ch > 0.3. Non-Fermi liquid temperature dependence 
is observed in chemical potential, specific heat, and uni- 
form spin susceptibility, at least in the region of inter- 
mediate temperatures 0.15 t < T < t, while at very low 
T < 0.15 t we cannot exclude a possible QP charac- 
ter of low-energy excitations, as signified by nearly con- 
stant Xo(T) and Cy oc T. We observe in addition a very 
non-monotonous Ch(fJ>,T) dependence, which we relate 
to a transition from a small hole-pocket-like to a large 
electron-like FS around doping c* h ~ 0.15. Entropy was 
found to be consistent in magnitude with experiments on 
cuprates, again showing maximum at Ch ~ c* h . Finally, 
the calculated Xo(r, Ch) is consistent with experiments, 
confirming the t — J model as a suitable one for the ex- 
planation of normal-state properties of cuprates. 
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FIG. 1. (a) Hole density c h as a function of T at several 
values of hole chemical potential fi h in steps of 0.1 t. Dashed 
lines denote fih < Mh • (b) ch as a function of fiu at various T. 

FIG. 2. Entropy density s — S/N vs. hole concentration 
Ch at several T. 

FIG. 3. Specific heat per site vs. T at several Ch- Ch — 
was obtained for N = 26. 



FIG. 4. Uniform susceptibility xo vs. T at several c h in 
steps of 0.05. Ch — was obtained for N = 26. 
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